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A renormalization group approach is used to show that a one dimensional system of bosons subject 
to a lattice quench exhibits a finite-time dynamical phase transition where an order parameter within 
a light-cone is zero up to a critical time and becomes non-zero after this time. Such a transition is also 
found for a simultaneous lattice and interaction quench where the effective scaling dimension of the 
lattice becomes time-dependent, crucially affecting the time-evolution of the system. Explicit results 
are presented for the time-evolution of the boson interaction parameter and the order parameter for 
the dynamical transition as well as for more general quenches. 
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A fundamental and challenging topic of research is to 
understand noncquilibrium strongly correlated systems 
in general, and how phase transitions occur in such sys- 
tems in particular. While the theory of equilibrium phase 
transitions is well developed, and relies heavily on the 
renormalization group, the development of an equally 
powerful approach to study nonequilibrium phase tran- 
sitions is still in its infancy. Moreover, in any progress 
on this topic, it has always appeared that nonequilib- 
rium phase transitions have one aspect in common with 
their equilibrium counterparts, both occur by adiabati- 
cally tuning some parameter of the system, in the absence 
or presence of an external drive, and strictly speaking oc- 
cur only in the limit of infinite time (steady-state) 

In contrast, here we study a completely different kind 
of a noncquilibrium phase transition, one that occurs as a 
function of time. Employing a time-dependent renormal- 
ization group (RG) approach we study quench dynam- 
ics of interacting one-dimensional (ID) bosons in a com- 
mensurate lattice. This system in equilibrium shows the 
Berezinskii-Kosterlitz-Thouless (BKT) transition sepa- 
rating a Mott insulating phase from a superfiuid phase 
(Fig. [T])jS For the nonequilibrium situation we explic- 
itly show the appearance of a dynamical phase transition 
where there is a finite critical time before which an order- 
parameter is zero, and after which it is non-zero. Such a 
behavior has no analog in equilibrium systems. 

A dynamical transition in time was recently identified 
in the exactly solvable transverse field Ising model where 
the Loschmidt echo was found to show non-analytic be- 
havior at a critical time, whereas the behavior of the 
order-parameter was analytic^ In contrast here we iden- 
tify a situation where the order-parameter itself can show 
non-analyticities as a function of time. In addition we 
generalize the study of dynamical transitions to mod- 
els that are not exactly solvable, and to low-dimensions 
where strong fluctuations negate a mean-field analysis; 11 

We identify the dynamical transition by studying an 
order-parameter A(r,T m ) which due to the quench de- 
pends both on position r and a time T m after the 
quench. The phase transition is associated with a non- 
analytic behavior as a function of time T m on the value 



of this order-parameter spatially averaged within a light 
cone. Our results hold relevance not only to exper- 
iments in cold-atomic gases where system parameters 
can be tuned rapidly in time^ but also to conven- 
tional solid state materials where time-evolution of an 
order-parameter may be probed with high precision us- 
ing ultra-fast pump-probe^ and angle resolved photoe- 
mission spectroscopy^ 

We model the ID Bose g Luttinger liquid. 
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where —d x 4>/ir represents the density of the Bose gas, II 
is the variable canonically conjugate to (j), Kq is the di- 
mensionless interaction parameter, and uq is the velocity 
of the sound modes. We assume that the bosons are ini- 
tially in the ground state of ffj. The system is driven 
out of equilibrium via an interaction quench at t=0 from 
A'o — > K , with a commensurate lattice V sg also switched 
on suddenly, at the same time as the quench. This trig- 
gers non-trivial time-evolution from t > due to a Hamil- 
tonian Hf = Hfo + V sg , where if/o = Hi(Ko — > K) and 
V S g = / dxcos(j(j)), with g > 0, and A = ^ a short- 
distance cut-off. We assume that the quench preserves 
Galilean invariance so that uK=uqKq, however this is 
not critical for either the approach or the result. While 
7=2 for bosons, we keep it general so that the results 
may be generalized to other ID systems. 

In the absence of the lattice, the system is ex- 
actly diagonalizable in terms of the density modes, 

H i=T, P7 io u o\p\v p Vp and Hf = Ep^o u \p\PpPp where /3, rj 
are related by a canonical transformation. This fact has 
been used to study the dynamics of a Luttinger liquid 
exactly, and has revealed interesting physics arising from 
a lack of thermalization in the systemJ^r— To study the 
system in the presence of the lattice employing RG, we 
write the Keldysh action representing the time-evolution 
from the initial pure state \</>i) (hence an initial den- 
sity matrix p — \4>i){4>i\) corresponding to the ground 
state of Hi, Z K = Tr[p(t)} = Tr [e-^l&X&le**'*] 
=f'D [4>ci, 4>g\ e l ( So+Ssg \ So is the quadratic part which 
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FIG. 1. The equilibrium BKT phase diagram. Arrows con- 
nect the Hamiltonians before (Hi ) and after (Hf) the quench. 
A dynamical phase transition is found for case (d). 



describes the nonequilibrium Luttinger liquid, which at 
a time t after the quench is?^ 

/oo />oo />t ft 

dxi / dx 2 / dh / tfta(^(l) 4>*(2)) 
-oo J-oo JO JO 

G A \1,2) \ U cl (2)\ 

G^(l,2) -[G^G K G A l ] {1,2)) \<f> 9 (2)J [Z) 

where l(2)=(a;i(2),ii(2)) and (/>cl,q= J + with — /+ rep- 
resenting fields that are time/ anti-time ordered on the 
Keldysh contour.— The lattice potential is given by 
Ss 9 =^ JZo dxi Jo dh [cos{ 7 <Ml)} " cos{ 7 + (l)}] 

We define an order-parameter A m = (e lm7 ^ c, ( a: ' t -'), with 
to < 1 which is zero in the gapless phase, and non-zero in 
the gapped phase, and a two-point correlation function 
C o6 (l,2)=(e i T'M 1 )e-n''M 2 )) where a,b=±. In order to 
understand the framework of the RG let us study C a b 
for .g=0 but for the nonequilibrium Luttinger liquid 
(K ? K). Denoting l(2)=R + (-)§,T ro + (-)§, C ab 
depends both on the time-difference r as well as the 
mean time T m after the quench, and is translationally 
invariant in spaced C a b depends on three exponents 



K r , 



eq- — ,^ n eq—8-K (l + fr) ,K tr =\K (l - ^ . 

For simplicity, consider C a b at equal-time (r=0) and 
unequal positions. Then for short-times T m A -C 1 
after the quench but long distances r u/A, C a b 
decays in position is a power-law with the exponent 
K neq + Ktr^Kv/A (i.e., C ab ~ r^ 2 *"/ 4 ). Hence 
the short time behavior is determined completely by 
the initial wave-function. In contrast, at long-times, 
T m A 3> 1, C ab decays as a power-law but with a new 
exponent K neq (i.e., C a b ~ r~ ne ' 1 ). K tr governs the 
transient behavior connecting these limits. Further, at 
long times after the quench T m A 3> 1, C a b also becomes 
translationally invariant in time (hence independent of 
T m ). 

The RG in equilibrium sums the leading logarithms. 
We use the same philosophy to employ RG to study 
dynamics. In particular at short times (but long dis- 



In Irl 



tances), the RG will resum the logarithms 
whereas at long times (T m A ^> 1), it will resum the loga- 



rithms K neq In yj (r ± r) 2 . Our approach generalizes the 
use of RG to study quench dynamics near classical criti- 
cal points^ 2 - to quantum systems. 

Derivation of RG equations: We split the field 
0o,A into slow (^ A _ dA ) and fast (</>A_ rfAjA ) compo- 
nents in momentum space <fi± = <fi± + cj)±, and inte- 
grate out the fast fields perturbatively in g. Following 
this we rescale the cut-off back to its original value and 
rescale position and time R,T m — > ■A(R,T m ), where 
A' = A — dA. Following this we write the action as 
SS^ + 5S£ + 5S< where is sim- 
ply the quadratic action corresponding to Hqj with the 
rescaled variables, Sf g is the rescaled action due to the 
lattice, while SSqt n are corrections to 0(g 2 ). 
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dR dT m [cos^<(R,T m ) 
) Jo 

- cos 70$ (R,T rn )} e-^^T) 

a 2 i 2 dA f°° f tA/V2 
SSf = Zf- -j- J dRj^ dT m [-I R {T m ) (d R i 

x (d R cf><) - I ti (T m ) (d Tm ^) (d Tm cf><)] 
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ig 2 j 2 dA 
"IT" T 



ff 2 7 2 dA 



(4) 
(5) 



3 ptK/s/2 

dR dT m (cp<) I Teff (T m ) 

a Jo 
oo ptA/V2 

dR / dT m <\>< [d Tm <j><] I n (T m ) (6) 

-oo JO 

Eq. © shows that the scaling dimension of the lat- 
tice depends on ([^(T m )] ), and in particular is time- 
dependent. To leading order (g—0) ^([^(^m)] 2 ) = 



A 



A 



ne ? 1 l+(2T m A) 2 



ss< 



shows that the quadratic 



part of the action acquires corrections which are also 
time-dependent~2, 5S^ T . indicates the generation of 
new terms such as a time-dependent dissipation (77) and 
a noise (rfT e ff) whose physical meaning is the genera- 
tion of inelastic scattering processes which will eventu- 
ally relax the distribution function££ These time depen- 
dent corrections lead to RG equations which depend on 
time T m after the quench. Defining t^ : = A a A =dln(/), 
I u ,k=Ir ± hr^i and dimensionless variables T m — > 
T m A,rj — > T]/A,T e ff — > T e ff/A the RG equations are, 
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= -T 



Ku d In I 
dr\ 



d(yT eff ) 

dlnl 



d\nl 
2r]T eff 



= V 



irg 2r y 2 K 



itg 1 K 



lT eff (T„ 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 
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Note that T m not only acts as an inverse cut-off in that 
modes of momenta A > 1/T m dominate the physics at 
a time T m , 24 i 25 but it also governs the crossover from 
a short time behavior where the physics is determined 
primarily by the initial state, and a long time behavior 
characterized by a new nonequilibrium fixed point. This 
crossover is most easily seen from Eq. ([7]) where the seal- 



ing dimension of the lattice e(T m ) 

depends on time as follows: at short times (T m <C 1) it 

2 Tv- 
is (— 2 + 7 4 ), and hence depends on the initial wave- 
function, at long times (T m 3> 1), a nonequilibrium scal- 
ing dimension (— 2 + K neq ) emerges. 

Above, lK,u,ri,T sJf reach steady state values at T m ^> 1, 
whereas for short times, they vanish as T m — > as ex- 
pected since the effect of the lattice potential vanishes at 
T m =0. For example, at short times I K ~ 0{T^)^L 
Eqns. ([5]) and © represent renormalization of the inter- 
action parameter and the velocity. The effects of the lat- 
ter being small will be neglected, and in what follows we 
set u=l. Eqns (fT0|) . (fTTj) show the generation of dissipa- 
tion and noise that represent inelastic scattering between 
bosonic modes££ 

In what follows we do an analysis for a time T m < \jr\ 
where 1 / 77 is the time in which the distribution func- 
tion first begins to change due to inelastic scattering. A 
perturbative calculation 5,6 shows that for small quenches 
(\K — K\ — > 0), and at steady-state, r\ ~ g 2 (K a — AT) 4 . 
Since ry <C 1 , one may easily be in the regime of T m 3> 1 
but T m r\ < 1 so that inelastic scattering may be ne- 
glected. At these times, and in what follows we will only 
use equations (J7J, flSJ) and (0). 

The behavior of the system is very different depending 
upon K ,K,g. We discuss four cases (see FigQJ. Case 
(a) is when the periodic potential is irrelevant at all times 
after the quench, case (b) is when the periodic potential 
is always relevant, case (c) is when the periodic potential 
is relevant at short times, and irrelevant at long times, 
while case (d) is when the periodic potential is irrelevant 
at short times and relevant at long times. For case (d) we 
show that an order-parameter behaves in a discontinuous 
way in time. There is a critical time T m before which 
the order-parameter is zero, and after which the order- 
parameter begins to grow indicating a dynamical phase 
transition. In contrast, for case (c), no discontinuous 
behavior of the order-parameter occurs. 

We use eo , <7o i ?m0 , Ao to denote bare physical values. 
We will present analytic solutions of the RG equations 
in the limits of short (T m oAo <C 1) times (but long dis- 
tances) and at long times (T m0 Ao ~> 1). For the for- 
mer the physics is primarily determined by the initial 
state. From Eq. we define an effective-interaction 

9eff { T m)=9\[^-^/ Ik {T m ) ^ J~j^ where g eff goes to 
zero as T m — > and reaches a steady state value for 
T m 1. Physically this implies that at short times 
the particles have not had sufficient time to interact, 
therefore however large g may be, any renormalization 
effects due to interactions is vanishingly small. The time- 



dependence of g e ff(T m ) and e(T m ) will be important for 
the results. 

Case (a), periodic potential always irrelevant 28 : 

This occurs for e(T m ) > and g e ff not too large (a 
condition to be made more precise when discussing case 
(d)). Here the periodic potential renormalizes to zero, 
and one recovers a gapless theory which eventually looks 
thermal at T m ^> l/r/£& The RG predicts how quantities 
renormalize in time and in particular shows that at long 
times the steady-state state is approached as a power-law 
with a non-universal exponent. 

Case (b), periodic potential always relevant^: 
This occurs for e(T m ) < 0. Thus we are always in the 
strong coupling regime. Here we integrate the RG equa- 
tions upto a scale l*(T m ) where the renormalized cou- 
pling is 0(1). Beyond this scale our RG equations are 
not valid, however the advantage of the bosonic theory is 
that at strong-coupling gcos(j4>) ~ g(l — r y 2 cj) 2 /2 + . . .) 
so that yfg may be identified with a gap. The physical 
gap/order-parameter is then given by A = ^fgjl* = l/l* . 
Since l*(T m ) depends on time, it tells us how the order- 
parameter evolves in time^ 

Let us first consider short times. Here we find that 
the order-parameter increases as A(X! m oAo <C 1) — 

(A T m0 .9o) L 4 ' ■ For a pure lattice quench (K =K) 
at the exactly solvable Luther-Emery point ^ we find 
A ~ 0(T^ Q ) in agreement with^°- The RG has allowed us 
to obtain results for the time evolution for more general 
quenches that take the system away from exact solubility. 

At long times after the quench, the scaling dimen- 
sion is 2 — K neq . Provided that AT m o 3> 1, we find 



the steady-state order-parameter, A ss — (<?e//,o) 2_JC ™ e ' ! ■ 
Compare this with the order-parameter in the ground 

state of Hf2- A eg =(<7 e //,o) 2 ^ K "i . Since K neq > K eq and 
9eff.o 1, the order-parameter at long times after the 
quench is always smaller than the order-parameter in 
equilibrium. The RG equations may also be solved at 
intermediate times 2 ^ t- T m o <C ir. Here we find, 



(A T„ 



2li£o . 



K„ 



?e//,0 



Thus at inter- 



mediate times the gap decreases with time if K neq > 

7 4 , or increases with time for the reverse case. Note 
that for Kq=K this intermediate time dynamics is ab- 
sent, instead the order-parameter reaches a steady-state 
at times T m oAo ^> 1. 

Case (c), periodic potential relevant at short 
times, and irrelevant at long times 28 . This oc- 
curs when e(T m ) changes sign from negative to posi- 
tive and g e ff is not too large. Here the short time 
behavior is the same as Case (b), however at long 
times, the order-parameter decreases with time as A ~ 

where A = 



( A T m0 ) 



- ' vv iiCic ^ . v5 2 //0 (oo)-eg(oo). 
Fig. [U summarizes the behavior of the order-parameter 
for cases (b), (c) and (d), the last case to be discussed 
next. The non-monotonic dependence of the order- 
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1/A„ X* 1/A - ^° 



FIG. 2. Time evolution of the order-parameter after the 
quench for cases (b), (c) and (d). Solid lines show short 
time behavior (T m o <C — ), an intermediate time asymptotics 

(xj ^ Tm0 ^ All) and a lon S time behavior T m0 » 
At short times for cases (b) and (c), the order-parameter in- 
creases as A k TmOi a t intermediate times it increases as T^ 
(decreases as T~q) when /C ne9 < ^ (_fr ne? > 2^a) for 
case (b) and eventually reaches a steady-state value A ss . For 
case (c) the order-parameter always decreases for T m o 2> j- 

as A oc T^q. For case (d) there is a non-analytic behav- 
ior at time T* n after which the order-parameter increases. 

9= ^ Ka , 0=0 _ Kncq \j ,j=l + M. Dashed lines 

are a guide to the eye for the crossover regimes. 



parameter in time is due to the time-dependence of 
the scaling dimension which physically leads to a sit- 
uation where quantum fluctuations are enhanced (sup- 
pressed) at a later time for e(T m = oo) > e(T m = 0) 
(e(T m = oo) < e(T m — 0)), causing the order-parameter 
to decrease (increase). 

Case (d), periodic potential irrelevant at short 
times and relevant at long times 28 : This occurs un- 
der two conditions. Either e(T m ) changes sign from posi- 
tive to negative during the time-evolution, or e(T m ) is al- 
ways positive, but g e ff(T m ) becomes sufficiently large at 
some time T^. The latter includes the case of a pure lat- 
tice quench (Kq=K). For either condition, at long times, 
the order-parameter reaches a steady state value, while 
at short times it is zero. This indicates a non-analytic 
behavior at a critical time T* n . 

Fig. [T] contrasts case (d) with the previous cases con- 
sidered where the order-parameter behaved analytically. 
The renormalized interaction parameter g e ff(T m ) is van- 
ishingly small right after the quench. For case (b), since 
infmitesimally small ff e // is a relevant perturbation, an 
order-parameter starts growing immediately after the 
quench. On the other hand for a quench corresponding 
to case (d), Fig.[T]shows that g e ff has to be larger than a 
critical value in order to be in the Mott-phase. Thus one 
has to wait some finite time before which renormaliza- 
tion effects become large enough for an order-parameter 
to grow. We now discuss this physics in a more quanti- 
tative manner, and for simplicity, consider only the case 



of the pure lattice quench. 

Let us suppose T m0 Ao S> 1. Here the RG equations 
are solved in two steps, one for 1 < I < T m oAo, and the 
other for T m oAo < I. For the first step, since Ik varies 
slowly at long times, eventually reaching a steady-state 
value, we may assume |||^| < |- Thus the RG 

equations are the conventional ones of the equilibrium 
BKT transition ^gf =-g effe , ■^ n =-g 2 eff . For the sec- 
ond step, (A T m0 < I), since Ik ~ T^, the RG equations 

2 _ 

become -^ff = -g eff e, ^ = where 1= 

The solution shows that there is a critical time such 
that g e ff is irrelevant before this time, and is a relevant 
perturbation after this time. We find, 

A T^= ei arcta„(ia)-i arctan % jfi = ^ _ £ 2 ^3) 

The deeper one quenches into the Mott-phase, the shorter 
is T* n . Moreover, T* n is longest along the critical 
line g e fffi=e$. By identifying a length-scale at which 
g e ff(l*) ~ 1, we find that the order-parameter grows as 

A~ 6(T m0 - T^)—^—[g eff {l = A T m0 )] ^14) 

where f 2 = | i , l , =Tm) 1 l -• 

An important question concerns the spatial variation 
of the order-parameter. Quenches in gapless systems are 
associated with light-cone dynamics where two points a 
position R apart get correlated after a time T m ~ i?£L 
For our case any two points separated by R > T m will 
behave primarily like the initial state with power-law cor- 
relations in position determined by Kq. The predictions 
for the order-parameter made above is for a region within 
a light cone R < T m . The dynamical transition at is 
associated with the appearance of order in regions of size 
R* ~ T* n , after which the ordered regions will begin to 
grow in size. 

In summary, employing RG we have identified a novel 
dynamical phase transition in a strongly correlated sys- 
tem where an order-parameter is zero until a critical time 
and increases after this time (Eq. (IHll ). The order pa- 
rameter shows rich dynamics both at the transition as 
well as for more general quenches (Fig. [2]). Identifying 
similar dynamical transitions in higher dimensions where 
thermal fluctuations are less effective in destroying order 
is an important direction of research. 
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Supplementary Material 



I. ELEMENTS OF THE QUADRATIC ACTION 
AND OUTLINE OF THE RG PROCEDURE 



The action for the non-equilibrium Luttinger liquid is 
S = I d Xl I dx 2 I dtt I dt 2 (<f>* cl (l) <t>* q {2)) 

(15) 



G A \1,2) 
G R \l,2) -[G-^GkG-/] (1,2) 

Denoting 1 = xi, ti; 2 = £2,^2 
= S(xi - x 2 )S(t 1 - t 2 ) 

while 



M2) 

M 2 ) 



irKu 



±1 



±iS ~ °x 1 \ 



dp 
— < 

P 



(16) 



cos(up(x - y)) cos(u\p\(t 1 - t 2 )) 



.Ka 



1 



ft' 2 
^0 



dp 



2 V *6/Jo P 
cos(up(x - y)) cos(u|p|(ti + t 2 )) 



where 



a 



(17) 



(18) 



is a short-distance cut-off. 

In doing the RG, the fields <fi are split into slow fields 
(</> < ) that have Fourier modes in momentum space be- 
tween 0, A — dA, and fast fields (</> > ) that have Fourier 
modes in momentum space between A — dA, A. Thus 
4> = 0< + > . The fast modes are integrated out per- 
turbatively in the periodic potential. In doing so we use 
the fact that the correlator for the field <fi (Go,a) is re- 
lated to the correlators for the slow dA ) and fast 
fields (Gl_ dA A ) as follows G ,a = G< A _ dA + G A _ dAA , 
so that— 



^A— <2A,A — ~dJ\~ 



(19) 



Following this we rescale the cut-off back to its origi- 
nal value and in the process rescale position and time to 
%i, U — > 4\(xi,ti) where A' = A — dA. 



II. EXPRESSION FOR C+_(r, T m , t) 

The two-point correlator defined as 

C + _(r,T m ,r) = ( e ^+(--^ ro +r/2) e -i0_(o,T m -r/ 2 ) )(2O) 



is for g — given by 
C+-(r,T m ,r) 



\/ o? + (ut + r) 2 \f a 2 + (ut — r) 2 



y/a 2 + {2u(T m + t/2)} 2 y/a 2 + {2u(T m - r/2)} 2 
y/a 2 + (2uT m + r) 2 y/a 2 + (2uT m - r) 2 

xe - 4 K og [ta„- 1 (i^±i:)+tan- 1 (^)] 



III. EXPRESSIONS FOR I R , I u , I n ,I Tet} , Ik, U 



-I Ktr 



(21) 



I R (T m A)-- 



f2T m A 

df / dff 

I -oo Jo 

Im [B(f,T m A,?)} 

/oo r2T m A 
df / dff 2 
-oo Jo 

Im [B(f,T m A,f)] 

/oo p2T m A 
df / df 
-00 J-2T m A 

Re[B(f,T m A,f)} 

/oo p2T m A 
df / dff 
-00 J-2T m A 

Im [B(r,T ro A,f)] 

/oo r2T m A 
df / df (f 2 + f 2 ) 

-00 JO 

Im[B(r-T m A,r)] 

/oo p2T m A 
df / df (f 2 - f 2 ) 

-00 JO 

Im [B(f,T m A,f)} 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



where Re[B] = (B + B*)/2, ha[B]=(B - B*)/(2i), and 
B(f, T m A, f ) = C+_ (f , T m A, f )F(r, T m A, f ) (28) 
with C_| (f , T m A, f) = 

^ e i* + (r,f+r m A/2) e -i0_(O,f-r m A/2^ < This quantity 

within leading order in perturbation theory is given by 
Eg. 1211 which we rewrite in dimensionless variables, 



C+_(f,T m A,f) = 



v/i + (f + f) 2 v/i + (f - ff 



y/l + {2(r m A + f /2)} 2 + {2(T m A - f/2jp 
v/l + (2T m A + f) 2 ^i + ^A-r) 2 

x g-iKeg [tan" 1 (f +r)+tan" 1 (r-f )] 



(29) 
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while F is given by 
F(r,T m A,f) = K neq 

■K. 



l + (f + f) 2 l + (f-r) 2 



1 + (2T m A + ff 1 + (2T ro A - f) 2 



-iif P 



1 + (f + ff 1 + (f - f) 2 



and 



K, 



j 2 K 



eq 



7 2 ( K 2 

K tr = —K 1 , 

8 V 



(30) 

(31) 
(32) 

(33) 



IV. SHORT TIME BEHAVIOR OF I K 

Ik depends on the two-point function C_| which may 

be computed explicitly within perturbation theory. We 
discuss its short-time behavior (T m A <C 1) here. In this 
case the integrand in Eq. [27] may be Taylor expanded in 
r so that 



I K (T m « 1) 



2K, 



dr 



2T„ 



drl r 2 



1 + r 2 



1 2 K Q 



(1 + r 2 ) 2 



eq 



1 2 K 



- 1 



0(t 3 



4 ^T 2 m + 0(T? n ) 



r(i 



(34) 

Note that Ik generically behaves as Ik ~ at short 



1 so 



times. However at the Luther-Emery point 
that the term of O(T^) vanishes, and the leading behav- 
ior is I K ~ O(T^). 

It is important to note that the Luther-Emery point is 
quite deep in the Mott phase, and is accessible perturba- 
tively in g only at short times. At long times, the slow 
power-law decay in time (and in space) leads to infrared 
singularities that makes perturbation theory break-down. 
In this paper we do not attempt to discuss the long-time 
behavior at the Luther-Emery point. 



V. SOLUTION OF THE RG EQUATIONS 

If we are interested in times that are such that T m < 
1 /r], where 1 jr\ is the time in which the distribution func- 



tion of the bosons changes considerably due to inelastic 
scattering, we may simplify the RG equations to 



dg 
dlnl 



= 9 



2 - 



K 

dK- 1 



K tr 



1 + 4T1 



din I 



1 2 
7T5 7 



dT m 
d\nl 



lK{T m ) 



= -T„ 



(35) 

(36) 
(37) 



We will use the notation go , T r „o, Ao as the bare physical 
values. Let us define the scaling dimension of the lattice 

as 



e(T m ) 



K 



K f , 



1 + 4T 2 



(38) 



We now consider four cases separately: Case (a) is when 
e(T m ) > at all times. Case (b) is when e(T m ) < at 
all times. Case (c) is when e(T m ) < at short times and 
e(T m ) > at long times. Case (d) corresponds to the 
dynamical transition and occurs for two possible cases, 
one where e(T m ) > at short times and e(T m ) < at 
long times. The second is when e(T m ) > at all times, 
but g e ff{Tm) becomes sufficiently large at long times. 

While one may always solve the equations [35l [36l [37] 
numerically, we can also obtain analytic solutions in two 
limits, one for short times T m oAo <C 1 and the other for 
long times T m oAo 3> 1. In these two limits a simpli- 
fication comes from the fact that 7if(T m ) for T m 1 
changes very slowly, as a power-law with T m , whereas 
at short times it changes as X^. For long times, g 
changes more rapidly with I than Ik does, so we will 
neglect the I dependence of Ik at long times. How- 
ever at short times, we will retain its I dependence via 
Ik (I™ « 1) - cT 2 A 2 /; 2 . 



A. Case (a) 

Case (a) is when e(T m ) > at all times, and g is not 
too large (this statement will be made more precise when 
discussing Case(d)). In this case, the periodic potential 
is always irrelevant. 

Let us first consider the short time solution. Here e = 
- 2. Further we write I K (T n 
the RG equations become 



« 1) 



dg 
din I 

dR^ 1 



7T3 2 7 2 / cT, 



2 A2 
raO A 



dlnl 



\ I 2 



cT^ . Thus 

(39) 
(40) 



Solving this upto I — ¥ oo (which is justified as g de- 
cays rapidly with I), the solution is g* — and 



K* 



_^ , , - £ + O{T 2 n0 k 2 ). Thus the renor- 

malized interaction parameter K* changes quadratically 
with time at short times from its bare value of K. 
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For long times T„ l0 Ao 3> 1, the RG equa- 
tions have to be integrated in two steps, one 
for T m 3> 1 where the RG equations become 
^=g (2 - K neq ) , J* (T m0 A /0. De- 

fine e = K neq — 2 where de = \j^-dK. We will now use 
the fact that the I dependence of /^(T m oAo/Z) is much 
weaker than that of g at long times i.e, ( 
1 9 ^ We also define an effective-interaction, 



1 1 din I K I 



9eff(Tm0^o) = 9 



^VWM^^ (4i) 



In terms of these new variables, in the long time limit, 
the RG equations are 



d 9eff 

d\nl 
de 



d\nl 



~ e 9eff 



-9eff 



(42) 
(43) 



The solution of the above equations are well known^ 

A 



e(lnZ) = 



tanh 



gQnl) 



A In I + tanh 
A 



-1 ( A 

to 



sinh 



A In 1 + tanh" 1 



where 



A = yJel(T m = oo) - g 2 fffi (T m = oo) 
eo(r m = oo) = if, 



neq 



(44) 
(45) 



(46) 
(47) 



We solve the above equations upto a scale Zi=T m oAo. 
Following this, for larger values of I, the RG equations 
change due to a change in the scaling dimension and also 
because Ik ~ T^. In order to smoothly connect the so- 
lutions of the RG equations at short and at long times 



we will make the ansatz, Ik 



/fr (exc- 



u- 



se) that in the second step, the RG equations change to 

dg.ff _ ( l 2 K a 2\ de _ 2 Tj A 2 , 

d\nl — \ 4 Z I 9eff, d \nl— 9 e ff P ' 1 n C SOIU 



tion of these equations with the initial conditions corre- 
sponding to Eqns. l44l l45l evaluated at 1% = T m oAo give, 



e*(A T m0 » 1) 



.4 



tanh 




Aln(A T m0 ) 
A 



tanh 



sinh 



Aln(A r m0 ) +tanh" 




Thus at long times a gapless phase is recovered, where 
the interaction parameter e* reaches its steady-state 
value of A as a power-law with a non- universal exponent, 
A T m0 >l i A , /n ( l 



-> A + O 



AnT„, 



2,1 



B. Case (b) 

This occurs for e(T m ) < 0. Thus we are always in 
the strong coupling regime. Here we integrate the RG 
equations upto a scale l*(T m ) where the renormalized 
coupling is 0(1). Beyond this scale our RG equations are 
not valid, however the advantage of the bosonic theory 
is that at strong-coupling gcos(-f<f)) ~ g(l — r y 2 (j) 2 /2 + 
. . .) so that y/g may be identified with a gap or order- 
parameter. The physical gap/order-parameter is then 
given by y/g/l* = 1/7*. Since l*(T m ) depends on time, it 
tells us how the order-parameter evolves in time. 

Let us first consider short times. Since g e ff(T m oAo <C 

1) ^ | 7 ?° — 2 1, we may approximate the RG equations 



de 
d\nl 







d 9eff ( I 2 Kg 



(49) 
(50) 



Integrating upto a scale I* such that g e ff=l, and using 
the short-time behavior of g e ff,o ~ goTmO, the effective 
gap is found to increase as 



A(T m0 A « 1) ~ (A T m0 g ) 
An order-parameter defined as 



is related to A = p- as follows, 



(51) 



(52) 



(53) 



A m (T m ) ~ [A(T m )] J 

where K(T m ) = AT V T m « 1 and K(T m ) = 
-^K neq V T m ^> 1. For a lattice quench (K =K ) at the 

exactly solvable Luther-Emery point (K eq — 7 4 = 1), 
Eq. El shows that I K ~ Tf n . Thus g eff ~ O(T^ ), so 
that A - O(T? n0 ) in agreement with RefM. 

At long times after the quench, the scaling dimension 
is determined by K neq . Here the RG equations are, 



de 







d 9eff 

dlnl 



n eq 



d\nl 

' 2 )9eff 



(54) 
(55) 



Integrating upto a scale where g e ff ~ 1, and provided 

that AT m0 ^> 1, we find the steady-state gap/order- 

i 

parameter, A ss — (g e ff,o) 2 ~ Kn '"' ■ Compare this with the 
gap in the ground state of A eq =(g e f ffl) 2 ~ K " q ■ Since 
K neq > K eq and g e ff,o "C 1, the gap at long times after 
the quench is always smaller than the gap in equilibrium. 

We now turn to intermediate times <C T m o <C 
here the RG involves two steps, the first in which the scal- 
ing dimension is determined by K neq and the integration 
stops at I* — AoT m o, and the second in which the scaling 
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dimension is determined by Kq, and the RG is termi- 
nated at g e ff = 1. Here we find (again taking ~ 0), 



A = 



(AoT m o) 



--K n 



lefffi 



Thus at inter- 



mediate times the order-parameter decreases with time 
if K neq > 7 ^ a , or increases with time for the reverse 
case. Note that for K =K this intermediate time dy- 
namics is absent, instead the order-parameter reaches a 
steady-state at times T m0 A ^> 1. 



C. Case (c) 



Case (c), periodic potential relevant at short times, 



rK 



< 



and irrelevant at long times. This occurs for 
2, K neq > 2 and g e f f not too large. Here the short time 
behavior is the same as Case (b). At long times, the RG 
proceeds in two steps. One where T m oAo > I > 1 where 
the solution is in equations|311l3S] For the second step the 



rK 



din l 



0. 



RG equations become -rrff = .9e//(2 
At the second step the RG equations are integrated upto 
a /* such that g e ff(l*) ~ 1. The solution of the order- 
parameter A ~ l/l* is found to be, 



A(A T m0 > 1) 
1 



AoT m o 



.4 



sinh 



Aln(A T m0 )+tanh- 1 



I 2 Kg 



(56) 



Thus at long times, the gap/order-parameter decreases 



to zero as A 



( A T m0 ) 



1+- 



D. Case(d) 

Here we will consider the case of the dynamical tran- 
sition where at short times the periodic potential is ir- 
relevant, whereas it becomes relevant at a finite non-zero 
time. Let us for simplicity consider this scenario for a 

2 K 

pure lattice quench so that K = K, with 7 -^— > 2 such 
that e > 0. This would naively imply that the periodic 
potential is irrelevant, however that is not the case when 
g e f f becomes large enough at some time . We demon- 
strate below how the dynamical transition occurs. 

Let us suppose we are at a time T m0 A ^> 1. At these 
times, the RG equations have to be solved in two steps. 
The first is when 1 < I < T m oA . Here, we may ne- 
glect the /-dependence of Ik by noting that it changes 
slowly in time at long times, in particular as a power- 
law, eventually reaching a steady-state. So assuming 
gl ^ffi ! < I TOT I in this regime, we get the follow- 



ing conventional RG equations 

d 9eff 

d\nl 
de 



d\nl 



9effC 



where 



9eff,o = 9o\/Ik(co)— I — 



(57) 

(58) 
(59) 



(60) 



The above equations may be solved by defining a con- 
stant of the flow 

D 2 =g 2 eff (l)-e 2 (l)=g 2 effi0 -e 2 (61) 



Integrating upto I = AoT m o we get, 

'e 



e(l — AoT m o)—D tan 



arctnu ( ) - Din (A T m0 ) 



g eff (l = A T m0 ) = ^D 2 + e 2 (l = A T mQ ) 



(62) 
(63) 



At the next step, the RG equations change because now 
Ik ~ l/l 2 . In order to smoothly connect the solutions 
of the RG equations at short and at long times we will 
make the ansatz, 



.'T m0 Ao , 
!k I — j — ) = Ik (oo) 



rp2 A 2 



The RG equations at the second step become, 
dg e ff 



d\nl 
de 
dhd 

I-- 



~9effe 

n 2 

9eff 



AoTtoo 



(64) 

(65) 

(66) 
(67) 



with the initial conditions at I = 1 given by Eqns [521 [Ml 
Equations [521 (Ml can be solved exactly. The solution is 



where 



e(l) = -1 + Vatan[ v / a(&-lnf)] 
9eff (l) = -y^sec [Vo(6-lni)] e lnI 



a = gl ff (l = l)-(e(l = l) + iy 
b = —= tan 1 



(68) 

(69) 
(70) 

(71) 
(72) 



The solution for g e jj{V) for different initial conditions are 
plotted in Fig. [3] There is a dynamical transition when 
g e /f (I = 1) = 9eff,c where 



leff,c 



- <U 2 (l = l) + 2e(l = 1) 



(73) 
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FIG. 3. RG flowof g e// (Z = 



for three different initial 



a re- 
conditions g e ff,i = g e ff(l = 1) and hence times, and for an 
initial = e(l = 1) = 0.1. Critical coupling or time is located 
at g eff>c = ^/2ei + e? = 0.458. 



The above critical coupling implies that there is a dy- 



namical transition at a critical time T* n such that 

g eff (l = AqO= ^(l = A T^) + 2e(Z = A T;„)(74) 



From Eq. 1631 the above implies 

D 2 = 2e(l = A T^) 
Using Eq. we solve for T^, 

AqT* = arctan (Tr)-lJ arctan § 



(75) 



(76) 



As before, we may identify the order-parameter/gap as 
the length scale I* at which g e ff{l*) ~ 1- This gives us 
the following result for how the order-parameter grows 
after the critical time, 



= 0(T m o ~ T* n )- 



1 



[ 5 e//G = A T m0 )]«^)-<^o) 

~ 6(T m0 - T*)— L- [g eff (l = A T mQ )}^d^j (77) 
where f 2 = l dei i =Tm u , • 

I ST^ l T m = T^\ 
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